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We describe generalizations of the Pauli group, the Clifford group and stabilizer states for qudits 
in a Hilbert space of arbitrary dimension d. We examine a link with modular arithmetic, which 
yields an efficient way of representing the Pauli group and the Clifford group with matrices over 
Zd. We further show how a Clifford operation can be efficiently decomposed into one and two-qudit 
operations. We also focus in detail on standard basis expansions of stabilizer states. 
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I. INTRODUCTION 



We study stabilizer states and Clifford operations for 
systems built from qudits (systems with a d-dimensional 
Hilbert space). We work in a matrix framework using 
modular arithmetic, generalizing results for qubits from 
Ref. Q . We put special emphasis on the less studied case 
where d is not prime. 

The stabilizer formalism has already proved to be use- 
ful in many applications such as quantum error correc- 
tion, entanglement distillation and quantum computa- 

tion iiai. 

The n-qudit generalized Pauli group and Clifford 
group and the related concepts of stabilizer codes and 
states have been studied in various levels of detail in a 
number of papers 0, 0, S, 0, E3, [HI E3, 113, Q, II^I 

Our motivation is not so much the study of stabilizer 
codes and their error correcting capacities, but the study 
of mathematically interesting states and operations that 
could play a role in quantum algorithms. Although it is 
well known that building quantum algorithms with sta- 
bilizer states and Clifford operations only is not sufficient 
to disallow efficient simulation on a classical computer, 
we think it is likely that the rich structure of this formal- 
ism will play a role in future quantum algorithms. Due 
to this focus, we pay attention to describing and realizing 
Clifford operations in more detail than is usually needed 
for coding applications. (To specify a Clifford operation 
"completely" (that is, up to only a global phase), one has 
to specify the image under conjugation of 2n independent 
Pauli operations including the resulting phase, whereas 
to realize an encoding operator for a fc-dimenional code, 
only k images are needed and the phases are of minor 
importance.) 

Next to presenting known results in an often different, 
and in our opinion practical language, we also present 
results not contained in the references above. 

We give a description of an rt-qudit Clifford operation 
by a 2n x 2n matrix C with entries in and a 2n- 
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dimensional vector h with entries in 1*2(1 and derive nec- 
essary and sufficient conditions for C and h to define a 
Clifford operation. We give formulas for multiplying and 
inverting Clifford operations represented in this way. 

We present a decomposition of a general Clifford oper- 
ation, specified in full detail by a matrix C and h, into a 
selected set of one and two-qudit operations, by thinking 
in terms of matrix manipulations on C and h. 

We also focus in detail on the standard basis expan- 
sion of stabilizer states. From Ref. 0, Q, formu- 
las can be derived describing the standard basis expan- 
sion of graph states by means of a quadratical form. In 
Refs. [lj, [l5j this is done for the case when the 1-qudit 
configuration space { 1 , . . . , d} is given the structure of 
a finite field. In Ref. this space can be any finite 
Abeli an g roup. In this paper we consider cyclic groups. 
Refs. [lj, [lj| state the equivalence of graph states and 
general stabilizer states. In Ref. 0] this equivalence is 
to be understood as local Clifford equivalence. That is, 
any n-qudit stabilizer state (with a field as 1-qudit con- 
figuration space) can be transformed into a graph state 
through the action of n one-qudit Clifford operations. In 
our setting however, as we are not focusing on codes, we 
want a description of the original stabilizer state (without 
the local Clifford operations) as well. In Ref. another 
notion of equivalence between graph states and stabilizer 
states is used (introducing the concept of auxiliary nodes 
in the graph). As a result the standard basis expansion 
of a general stabilizer state is not described directly but 
as a sum of a large number of states. Moreover, for the 
case where the configuration space is not a field (in our 
case that is when d is not prime) not all stabilizer states 
are equivalent to graph states but an extra condition has 
to be imposed. In the present paper we work with a more 
general description of stabilizer states without this extra 
condition (described below by matrices S with possibly 
more than n columns) and we give a direct description 
(without sum) of the standard basis expansion of general 
stabilizer states. We believe that standard basis expan- 
sions of stabilizer states can be an essential ingredient in 
understanding the action of non-Clifford operations on 
stabilizer states. 

This paper is structured as follows. Definitions of gen- 
eralizations of the Pauli group and the Clifford group for 
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qudits are given in section [HJ together with their ma- 
trix representation. Special Clifford operations, that are 
of particular interest in the decomposition of a Clifford 
operation, are discussed in section ITTT1 An efficient de- 
composition of a Clifford operation on n qudits into a 
selected set of one and two-qubit Clifford operations, is 
explained in section llVl In section[V] we define stabilizer 
states of n qudits and show the expansion in the stan- 
dard basis can be described with linear and quadratic 
operations. 

In the following, by A = B mod d we mean that all 
corresponding entries of matrices A and B are equal mod- 
ulo d, where d is an integer different from 0. We will also 
write a — b mod c with a, b and c vectors, as a shorthand 
notation for a t — bi mod Cj, for every i = 1 . . . ri- 



ll. THE GENERALIZED PAULI GROUP AND 
CLIFFORD GROUP 

In this section, we discuss the description of the gener- 
alized Pauli group on n qudits and the generalized Clif- 
ford group in modular arithmetic. Generalizations of the 
Pauli group to systems of arbitrary dimensions are dis- 
cussed in Refs. jatllal- The Clifford group is defined as 
the group containing all unitary operations that map the 
Pauli group to itself under conjugation. 



A. The generalized Pauli group 

Let d be the Hilbert space dimension of one qudit. We 
define unitary operations and as follows 



where U 



Addition in the argument of XZ 



Z {d) \j) = <J\$), 



(1) 



where j £ and uj is a primitive <i-th root of unity. Ad- 
dition in the ket is carried out modulo d. Tensor prod- 
ucts of these operations will be denoted as follows: for 

v 

D.aieZJ and a := 

a w 



£ 1r d n , we denote 



XZ(a) :=X Vl Z Wl ® ...®X Vn Z Wn . (2) 
From and it follows that, for x £ 1%, 

XZ(a)\x) =lo wTx \x + v). (3) 

We define the Pauli group V„ on n qudits to contain all 
d 2n tensor products with an additional complex phase 
factor ( 5 , where £ is a square root of cu and 6 £ Z 2 d- In 
the following, we will omit the superscript (d) and refer 
to the generalized Pauli group simply as Pauli group. 

Multiplication of two Pauli group elements can be 
translated into operations on vectors in Z^ n as follows: 



On On 
In n 

is done modulo d, and addition in the exponent of ( is 
done modulo 2d. Eq. (0J yields the commutation relation: 



XZ(a)XZ{b) 



Pb 



XZ(b)XZ(a) 



where P — U — U mod d. 



(5) 



(6) 



Note that the order of XZ(a) divides d unless d is an even 
number and a T Ua is odd. In the latter case the order is 
2d. Indeed, with one can easily verify that XZ{a) d = 
(■d(d-i)a Uaj^ introduction of a phase <^ 5 rather than 
uj s is only necessary when d is even. Simplifications for 
odd d are considered in Appendix ICl 



B. The generalized Clifford group 

We now define a generalization of the Clifford group 
on n qudits in an analogous way as for qubits. A Clifford 
operation Q is a unitary operation that maps the Pauli 
group on n qudits to itself under conjugation, or 

Because QXZ(a)XZ(b)Q^ = (QXZ(a)Q^)(QXZ(b)Q i <), it 
is sufficient to know the image of a generating set of the 
Pauli group in order to know the image of all Pauli group 
elements. Q is then defined up to a global phase factor. 
This can be seen as follows. Suppose that two Clifford op- 
erations Qi and Q2 give rise to the same image for every 
Pauli group element, or: for every A £ V n ■ QiAQ\ = 
Q 2 AQ\. It follows for every A that Q\QiA = AQ\Q x . 
The only unitary operations that commute with every 
single Pauli group element are multiples of the identity 
jlq . which completes the proof. We take the generat- 
ing set of the Pauli group to be XZ(Ek), k = 1, . . . , 2n, 
where E/~ are the standard basis vectors of . We de- 
note their images under conjugation by Q as ( hk XZ(Ck). 
We will assemble the vectors Ck as the columns of a ma- 
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( 5 XZ{a)CXZ(b) = C 



S+e+2a T Ub 



XZ(a + b), (4) 



trix C £ and the scalars hj~ in a vector h £ Z 2d 

The image ( € XZ(b) of Q & XZ(a) under conjugation by Q, 
where a is an arbitrary vector in Z^™, can be found by 
repeated application of QJ. This yields 

b = Ca mod d, 

e = 8+ (h-V diag {C T UC)) T a+ 

a T (2V upps (C T UC) + V diag {C T UC))a mod 2d, 

(7) 

where Vdiag(Af) is defined as the vector containing the 
diagonal of M, PdiagC-^) the diagonal matrix with the 
diagonal of M and V U pps(M) the strictly upper trian- 
gular part of M. The Clifford operation Q is (up to a 
global phase factor) completely defined by C and h. Note 
that the rhs of J7J is calculated modulo 2d, although it 
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contains matrices over Z^. It can be verified that every 
entry modulo d in the expression is multiplied by an even 
factor. 

We can compose two Clifford operations Q and Q 1 , 
which again yields a Clifford operation Q" = Q'Q. To 
find its corresponding C" and h" we have to find the 
images under the second operation of the images under 
the first operation of the standard basis vectors. By using 
Q , we get 

C" = C'C modd, 

h" = h + C T ti + V dllig (c T (2V upps (C ,T UC')+ 

v diag {c> T uc>))c)~c T v dmg (c< T uc') (8) 

mod 2d. 

The inverse of a Clifford operation Q defined by C 
and h is defined by C and h', where 

C = C'- 1 mod d, 

h' = -C- T [h + V d ^(c T {2V^ ps {C- T UC- l ) + 

v^, s {c~ T uc- l ))c) - c^wc^t/c- 1 )^) 

mod 2d, 



(9) 



which can be verified with (jSJ. M T is short for (M 1 ) ' 
We will show below that C" 1 = -PC T P mod d. 



C. Conditions on C and h 

Not all C £ z 2 d nx2n and h G Z;f d l define a Clifford 
operation. To see this, consider a Clifford operation Q 
with corresponding C and h. From the commutation 
relation (jSJ) it follows that C is a symplectic matrix, i.e. 
C satisfies C 7 PC = P mod d. Indeed, we have 

XZ{a)XZ(b) = oj aTpb XZ(b)XZ{a) 

QXZ(a)Q?QXZ(b)Q? = uj a ' ' pb QXZ(b)($QXZ(a)($ 

XZ{Ca)XZ{Cb) = uj aTpb XZ(Cb)XZ(Ca), 

where we omitted global phase factors on the lhs and rhs, 
as they cancel each other. Also, 



XZ(Ca)XZ(Cb) = lu 



a 1 C 1 PCb 



XZ(Cb)XZ{Ca). 



Since this holds for every value of a and b, it follows that 
C is symplectic. Note that the inverse of a symplectic 
matrix C is simply C" 1 = —PC T P mod d. Secondly, h 
satisfies 

(d-l)V diag (C T UC) + h = mod 2, (fO) 
for C hk XZ{C k ) = QXZ{E k )Q^ has, like XZ(E k ), order d. 

With Q we have (C hk XZ{C k )) d = ^((d-l^^+A*)^ 
and it follows that l|f Ofl is satisfied. We will prove below 
that every symplectic C and h satisfying (|f 0|) define a 
Clifford operation Q. 



III. SPECIAL CLIFFORD OPERATIONS 

In this section we present a number of special Clifford 
operations and their defining C and h. These will be 
of particular interest for the decomposition of an arbi- 
trary Clifford operation into one and two-qubit Clifford 
operations. 

• The Pauli group elements XZ(a) are a special class 
of the Clifford operations. Note that, like for any 
Clifford operation, the global phase factor of a Pauli 
group element cannot be represented. Considering 
the images of XZ(E k ), it can be easily verified that 
XZ(a) is defined by 

C = I mod d 
h = -2Pa mod 2d. 

• A Clifford operation acting on a subset a C 
{1, ...,n} of n qudits gives rise to a symplectic 
matrix on the rows and columns with indices in 
a U (a + n), embedded in an identity matrix (that 
is, C kk = 1 mod d, for every k ^ a U (a + n) and 
C k i — mod d if k ^ / and k or I ^aU(aln)). 
Also h k = mod 2d if k £ a U (a + n). 

• Any invertible linear transformation of the config- 
uration space \x) — > \Tx) can be realized by a Clif- 
ford operation, with x £ Z^ and T e ZJ X " an in- 
vertible matrix modulo d. This operation is defined 
by 



C = 
h = 



T 

T- T 




mod d 
mod 2d. 



This can be verified by looking at the image of 

v 
w 



XZ(a), with an arbitrary a 



QAZ(a)Qt 



£ \Tx){x\\ "£u wT «\y + v)(y\ 



E u wTy \Ty + Tv)(Ty\ 

E ^ TT " ly iy+^)(yi 

Tv 
T~ T w 



XZ( 



)• 



As C T UC = U mod d, we see with Q that h = 
mod 2<i. Special cases of this class of Clifford oper 

ations are qudit permutations, with C = |J 



4 



where II is a permutation matrix, and the two-qudit 
SUM gate \x)\y) — ► \x)\x + y) with x,y G Z d , with 



C = 



10 
110 
1-1 
1 



3d d. 



Note that this operation is a natural generalization 
of the two-qubit CNOT gate. 

• The d- dimensional discrete Fourier transform 
\x) — > 4= J2k=o ^ kx \k) on one qudit, with x G Z<j, 
-1" 



is defined by C = 



1 



mod d and h = 



mod 2<i. We verify this in the same way as for the 
invertible configuration space transformation, now 

v 



with a = 



G Zj: QXZ{a)Qi 



U d E «-'i*x*i) 
- 1 E 



w^ly + «)<w| 



t(l/+D)+«)j/-xy 



t,y,xeZd 



= 1 e w (tfu,_B) * E ^i*)( a 

= u {x ~ w)v \x-w){x\ 

)■ 



w~ utu JiZ( 



As C T UC = -U T mod d, we see with Q that ft = 
mod 2d. This operation is the qudit equivalent 
of the Hadamard gate on one qubit. 

• Analogous to the qubit phase gate, a phase gate on 
one qudit can be defined as \x) — > C (x+d) \x), with 

x G Z,j- This operation corresponds to C = J ^ 

mod 2d. Indeed, for all 



mod d and h 



d + 1 




G Z2 : QXZ(a)Q^ 

£2ioy+(iH-u)(i/+t;+<i)-y(y+<9|^ _|_ 

^ E ^ (t,+, " )s iy+«>(yi- 



As C T £/C = 



1 
1 



mod d, v(v + d) must be 



equal to h 



- v 2 mod 2d according to 



J7J), which is the case for the given h 



IV. DECOMPOSITION OF A CLIFFORD 
OPERATION IN ONE AND TWO-QUDIT 
OPERATIONS 

In order to prove that any symplectic matrix C and 
h satisfying i|10fl define a Clifford operation, we will ex- 
pand an arbitrary symplectic C into symplectic elemen- 
tary row operations that can be realized as Clifford oper- 
ations on maximally two qudits at the same time. What 
is more, this decomposition is a worthy candidate as a 
practical realization of a Clifford operation. The possi- 
bility of this kind of decomposition into a selected set 
of one and two-qudit operations is briefly discussed in 
Ref. Q . Our scheme is related to the method of Ref. [f| 
in which Euclid's algorithm is incorporated in order to 
generate any one-qudit Clifford operation. 

First, we mention that the main problem is realizing 
C, not h, for once a Clifford operation Q defined by C 
and h is realized, we can realize Q' defined by C and an 
arbitrary h! satisfying IjlOl) by doing an extra operation 



(CP^) 



on the left or XZ [P 



yh' — h 



on the right of 



2 j u .i »uv i-i. ^ ^ 2 

Q. Note that as both h and h' satisfy l)10|l. h! — his even. 

We first give an overview of the elementary row opera- 
tions that we will use to transform an arbitrary symplec- 
tic matrix C into the 2n x 2n identity matrix /. As I is 
formed by left multiplication of such elementary row op- 
erations on C, a decomposition of C then consists of the 
inverses of these operations in reverse order. Since these 
operations act on maximally two qudits at the same time, 
they are defined by a symplectic 4 x 4 or 2 x 2-matrix 
embedded in the identity matrix as explained in the pre- 
ceding section. In the following, we will only show this 
part of the operations. 

Firstly, we consider some configuration space transfor- 

T 
T- T 

combine only rows from the same block (we call rows 
1 . . . n the upper block and rows n + 1 . . . 2n the lower 
block) and have a similar action in both blocks at the 
same time. For instance, we can switch two rows i and 
j in the upper block: at the same time, rows n + i and 
n+j in the lower block are also switched. This operation 
is defined by 



mations (of the form C 



). These operations 



C = 



10 
10 
1 
10 



3d d. 



Multiplying a row i with an invertible number r G Z<j 
results in multiplying the corresponding row n + i in the 
other block by r _1 . A number r G Z^ has an inverse if r 
and d are coprime, i.e. gcd(r, d) = 1. This operation is 
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r 

r- 1 

transformation we consider 



defined by C = 



The last configuration space 



is adding one row i multi- 
plied by an arbitrary factor g G to another row j. At 
the same time, row n + j, multiplied by —g, is added to 
row n + i. This operation is defined by 



C 



10 

g i o o 

f -g 
f 



mod d. 



Secondly, we will also need operations that combine 
rows of different blocks. Switching two rows i and n + i 
can be carried out by the discrete Fourier transform. Re- 
call that this operation is defined by C " 



f 



Af- 



ter switching of the rows, row i is multiplied by — f . By 
applying the inverse of the discrete Fourier transform, 
row n + i instead of row i is multiplied by — 1. Apply- 
m S J2 x tEZ d C 9xi ~ x+d ^\ x )( x \ (which is the same as applying 
the phase gate g times) on the i-th qudit, with g e Z^, 
results in the addition of row i multiplied by g to its 
corresponding row n + i, according to 



C 



1 
9 1 



mod d. 



We could introduce more row operations that define 
one or two-qudit Clifford operations, but the ones de- 
scribed so far suffice. Next we give a constructive way of 
transforming C into the identity matrix /. If we are able 
to transform C into C by transforming columns C± and 
C n +i into the corresponding columns E\ and E n+ i of I, 
it follows from the symplecticity of C that the first and 
n + f-th row of C are equal to the corresponding rows 
of I. We then have 



c 



o ... o 



C (U) 



... 



C, 



(21) 



... 



(-11 

°(12) 



... 



(-11 

°(22) 



Leaving the first qudit out, we can continue by transform- 
ing the second and n + 2-th column of C into the cor- 
responding columns of I, and so on. This recursive pro- 
cedure eventually leads to I. Now we only have to show 
how columns C\ and C n +\ are transformed into E\ and 
E n+ i. Let us first consider the case where the upper left 
entry C\\ has an inverse in Z^. Multiplying the first row 
by C-j^ 1 changes this entry to 1 . Next we add the first row, 
multiplied by — Cki, to row k, and this for k — 2 . . . n, 
setting the k-th entry of C\ to 0. The first column now 



has the form [10 • ■ ■ | C' n+11 C n+2 ,i • • ■ C 2n s] T ■ 
Now we add the first row multiplied by — C' n , j 1 to row 
n + 1, setting the n + 1-th entry of C\ to 0. The dis- 
crete Fourier transform on the first qudit changes C\ into 
[0 • • • | 1 C„+2,i • • • C2n,i] T - In the same way as for 
the upper half of C\ , we make zeros below the n + 1-th 
position. Note that nothing happens to the upper half, 
for all entries there are 0. Switching the first (now we 
use the inverse of the discrete Fourier transform) with 
the n + 1-th row again yields E\. We call the matrix 
made so far C" . From the symplecticity of C" it fol- 
lows that C^'+i.n+i = 1 mod d, and we can repeat for 
the n + 1-th column the same procedure we did for the 
first column. Note that none of the operations yielding 
E n+ i out of C" +1 will affect C'{ — E\, except the discrete 
Fourier transform and its inverse on the first qudit, but 
they cancel each other. Since the number of elementary 
operations for one column is 0(n), the total number of 
operations transforming C into / is 0(n 2 ). 

If the entry C\\ has no inverse modulo d, but there is a 
Cfei in the first row that does have an inverse, this entry 
can be switched into the first position by a permutation 
of two qudits and possibly the discrete Fourier transform 
on the first qudit. Note that it is possible that none of the 
entries of C\ has an inverse. Indeed, since C is invertiblc, 
the only restriction on one single column of C is that the 
greatest common divisor of all its entries has an inverse. 
For every two entries Cn and C n +i t \ or Cn and Cji from 
the same block, the gcd of these two can be formed in one 
of the two entries by recursively substracting a multiple 
of one row from the other following Euclid's algorithm 
[l6j . The other entry can then be made since it is a 
multiple of the gcd. A worst case scenario would be that 
all 2n combinations of 2n — 1 entries have a gcd that is 
not invertible. The procedure goes as follows 



gcd(Cii,C n+ i,i) 
gcd(C 2 i, C„ +2 ,i) 

gcd(C„i,C 2 „,i) 







C21 












Cn+2,1 




- Cln,\ - 





gcd(C 



ii- 



. C 2n> i) 



In this way, C is decomposed into 0(n 2 log(d)) elemen- 
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tary operations, as the computational complexity for 
finding the gcd of two positive integers less than d with 
Euclid's algorithm is 0(log(d)) 16]. 



V. STABILIZER STATES 

In this section we define stabilizer states for qudits of 
arbitrary dimensions. A stabilizer state is a state of an 
77-qudit system that is a simultaneous eigenvector, with 
eigenvalues 1, of a subgroup of d n commuting elements 
of the Pauli group, which is called the stabilizer S of 
the stabilizer state. The stabilizer state is completely 
determined by a generating set for S. The description 
of such a generating set in modular arithmetic provides 
an efficient tool of describing the stabilizer state and its 
behavior under the action of a Clifford operation. Finally, 
we give an expansion of an arbitrary stabilizer state in 
the standard basis. 



A. Definition and description in modular 
arithmetic 

A stabilizer state l^) is the simultaneous eigenvector, 
with eigenvalues 1, of a subgroup of d n commuting el- 
ements of the Pauli group which does not contain mul- 
tiples of the identity other than the identity itself. We 
call this subgroup the stabilizer S of \ip). A generating 
set for S consists of elements XZ(Sk), k = 1 . . .m, 
where Sk € and fk € lid- We will assemble the 
vectors Sk as the columns of a matrix S € Z^ nxm and 
the scalars fk in a vector / 6 ZJ^j. We call S a gener- 
ator matrix and / the corresponding phase vector that 
together define S. The fact that the elements of S com- 
mute is reflected by S T PS = mod d. We choose m 
to be the minimal cardinality of a generating set of S. 
Note that, as opposed to the situation for qubits, m can 
be larger than n. It can be verified that if m > n, the 
imposed condition in Ref. for a stabilizer state to be 
equivalent to a graph state, is not fulfilled. If d has only 
single prime factors, then m = n. If d has multiple prime 
factors, then n < m < 2n. A simple example for d = 4 
and n = 1 is the state 1/V2(|0) + |2)) with stabilizer 
{/, X 2 , Z 2 ,X 2 Z 2 }: in this case m = 2. We will describe 
below how to construct such a minimal generating set. 
The fact that S does not contain multiples of the iden- 
tity other than the identity itself implies that the phase 
vector / satisfies: 

Vr e Zf | Sr = mod d : (f - V diag (S T US)) T r+ 
r T (V diag (S T US) + 2V upps (S T US))r = mod 2d. 

(11) 

The description of S by S and / is not unique, as they 
represent a generating set for S. By applying an invert- 
ible linear transformation R £ Z™ xm to the right on S 
and transforming / appropriately, another generating set 



Q^XZ{S' k ) is formed. By repeated application of one 
finds 

S' = SR modd, 

/' = R T (f-V diag (S T US)) + 

V dmg (R T (2V upps (S T US)+ ( 12 ) 

V dilie (S T US))R}j mod 2d. 

We will refer to this as a stabilizer generator matrix 
change. 

If \ip) is operated on by a Clifford operation Q, defined 
by C and h, then Q\ip) is a new stabilizer state whose 
stabilizer is given by QSQ^ . By application of J7J, we 
can calculate an S' and /' for this stabilizer, resulting in 

S' = CS mod d, 

f = f + S T (h~V diag (C T UC)) + 

V di . e (s T (2V upps (C T UC)+ ( 13 ) 

V diag {C T UC))S^j mod 2d. 

We can construct a minimal generating set 
^ k XZ(Sk), k — l...m, for an arbitrary stabilizer 
S, given a generating set XZ{S[), I = 1 . . . m', for S 
using the Smith normal form (see Appendix ^ . This 
can be done as follows. The S[ are assembled in the 
matrix S'. Now we compute the Smith normal form 
F = KS'L of S', with K e Z 2 d nx2n and L € Z™' xm ' 
invertible matrices. S'L is just another generator matrix 
of the stabilizer. From the definition of the Smith normal 
form it follows that S'L is a generator matrix having 
a minimal number of nonzero columns. The rightmost 
to — m' columns of S'L that are zero (as S'L = K~ X F) 
can be omitted. We call this new generator matrix S 
and / is formed out of /' with (|12() . Note that no linear 
combination of the columns Sk of S is zero unless the 
coefficients in this linear combination are a multiple of 
the order of the columns, or, for k = 1 . . . m, 

if ^""^ rkSk = mod d, then rkSk = mod d. (14) 

k 

With this, the stabilizer phase condition (|llfl can be sim- 
plified to, for k = 1 . . . m: 

W k e Z d I r k S k = mod d : 
(r k - l)r k SlUS k + r k f k = mod 2d. 



B. Description of a stabilizer state with linear and 
quadratic forms 

We provide an expansion of an arbitrary stabilizer 
state in the standard basis for an n-qudit state. This 
is stated in the following theorem 
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Theorem 1 (i) If S G z 2 d nxm and f G Z™ de/me a 
stabilizer state \tp) as described above, then S and f can 
be transformed by an configuration space transformation 
\x) — > \T~ 1 x), with T G Z^ x ™, and a stabilizer generator 
matrix change R G Z™ xm info the form S' and f , with 



S' 



(16) 

where Q is a pseudo- diagonal matrix in Smith normal 
form and Q T B mod d is symmetric. Q and B are the 
left square n x n parts of Q and B. 

(ii) The state \ip) can be expanded in the standard basis 
(up to a normalization factor) as 



T -i o 
T T 


SR = 


Q 
B B 




Q 

B 


mod d, 




f' T = 


r r 


T ' 


mod 2d, 



H') = E c 



t T Mt+p T t 



\T(Qt + x*)) 



where M := QB mod d, 

V ■•= f ~ Vdiag(M) + 2B T x* mod 2d. 

If we define the n-vector q with entries qk 

\ d if Qkk — mod d 
\ Qkk if Qkk 7^ mod d 



k = 1 



. n, anc 



the 




Then x* G G g 



aq n is defined as the unique solution of 



B x = y mod q, 



(18) 



where y G 

f (d-q k )B kk +f' k 

\ 2 -4 



has entries 



Ilk 



mod qk, for k = 1 . . . n 
mod qk, for k = n + 1 . 



Note that from the stabilizer phase condition (|15Jl (choose 
r k := d), it follows that the numerators in the expressions 
for yk are even. An efficient way of solving (|18ll can be 
found in Appendix iBl 

A definition of the Smith normal form of a matrix G 
z «xm ig giyen - n Appendix El 

Proof: 

(i) We assume that S already has a minimal number of 



columns m as described above and we write 5* as 



(2) 



with 5(2) G Z' d ixm . Then we define Q as the Smith 
normal form of 5(1) with invertible transformation matri- 
ces T" 1 and R, i.e. Q = T^S^R. With B = T T S [2) R, 
this yields the expression for S' in i|16|) . According to 
(|12(l and 1)13(1 . / is transformed to /', yielding 

/' = i? T (/-V diag (5 T f/5))+V diag (i? T (27' U pp S (5 T C/5) 
+V diag (S T US))R) mod 2d. 



Note that 



T- 1 
T T 



U 



T- 1 
T T 



= U mod d. It 



follows directly from S T PS — mod d that Q T B is 
symmetric modulo d. 

(ii) We show that (|17fl is a simultaneous eigenvector 
with eigenvalue 1 of £^ fe XZ(S k ), k — l...m. Equiva- 
lcntly, the state 



W) ■= E c 



t J Mt+p J t 



|Qf + x* 



(19) 



is a simultaneous eigenvector with eigenvalue 1 of 
tf'*XZ(S' k ), k = l...m. First, note that in jJHJ), different 
values of f may yield the same basis state \Qt + x*), since 
Qt + x* mod d is periodic. The coefficient of \Qt + x*) 
in (|19fl displays the same periodic behavior: if Qt — Qt' 
mod d then f T Mf + p T t = t' T Mt' + p T t' mod 2d. It is 



(17) sufficient to check this for t' = t 



E k , k = l..,n, 



where E k are the standard basis vectors of Wi. We have 



(f + f k E k ) T M{t + ±E k ) + p T (t + ±E k ) - t T Mt ~ p T t 



<ii. 



ik ' 



= i((d- q k )B kk + f' k + 2Blx*) = mod 2d, 



Qk 

for k 

QT X * — (d~q k )B kk +f' k 



l...n. Indeed, from the definition of x*: 

k Ju — 2 mod q k , k = 1 . . . n, it follows 

that 2f k B T k x* = -£((d - q k )B kk + f' k ) mod 2d, k = 
1 . . . n. We made use of the fact that M = QB mod d 
M = QB + D mod 2d, where every entry of D mod 2d 
can be either d or 0, i.e. 2D = mod 2d. 

Next, we check for k = 1 ... n that l|19f) is an eigenvec- 
tor of C,^ k XZ{S' k ) with eigenvalue 1. We have 



Qk 

B k 



W) 



= Ec 



t T Mt+p T t+f k +2B k (Qt+x* 



\Qt + X* +Q k 



= £(t-E k ) T M(t-E k )+p T ' (t~E k )+f' k +2Bl(Q(t-E k )+x*) 

\Qt + x*) 

= Y,c TMt+pTt \Qt+x*) = w). 

Finally, C f ' k XZ(S' k ) acting on the left of ||TSJ| yields, for 
fc = n + 1 . . . m, 







1^) 



^( fTMf+pTf |Qf + x*) = |V')- 
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In Appendix [B] we prove that eq. (|18fl has a unique solu- 
tion x* € Gg. □ 

It is possible to remove all identical terms in the sum- 
mation of expression (|T7I) as follows. We define r as the 
number of nonzero diagonal elements of Q. We denote 
the upper left r x r-part of a matrix A as , the upper 
r-part of a vector a as a( r ) and the part of a below ai r ) 
as a( r ). Then l|17|l is equivalent to 



n; 



E 

teG, 



^t 1 M (r) t+p{ r) t I j, 



kr)t + X* [r) 
X U 



x Z d . Note that the normal- 



where G* := Z_d_ x 

izing factor is just the inverse of the square root of the 
number of terms in the summation, as each basis state 
is orthogonal to the others and occurs only once. Fi- 
nally, it is interesting to mention that, for an arbitray 5 
and / defining a stabilizer state \tp), we have (up to a 
normalization factor) 



= E c tTMt+pTt \s {1) t+ x '), 



i.ez 



where 5 



Sn 



(i) 

(2) 



M 



mod d, 



P 



f 



Vdiag(A/) + 25^)a;' mod 2d and x' = Tx* mod d, where 
T and x* are the same as in (|17f) . Yet, this formula has 
two disadvantages: first, to find x' , we still have to cal- 
culate the Smith normal form of 5m and second, in ltT7|) 
it is clearer which basis states have nonzero coefficients. 



!%* m such that 



F = KAL 



mod d 







with each fi a nonzero and with fi\fi + \ for 1 < i < 
r — 1. The fi are unique up to units. Uniqueness of F 
can be ensured by specifying that each fi should be a 
positive divisor of d in Z. There exist fast algorithms for 
computing the Smith normal form |17|. 



APPENDIX B: A UNIQUE SOLUTION OF JT51) 

Here we prove that eq. <|18|) B T x = y mod q has a 

. x Z gn . We rewrite 



unique solution x* G Gq := Z 9l x 
(|18ll as the following system of equations: 



7 , -^i?" 
i=l 



% mod q.,-, j = 1 



(Bl) 



If, for fixed j, the Bij and have a common factor, then 
also ijj must be a multiple of this factor, otherwise there 
9o : = gcd(Sy , . . . , Bnj^j) and 



is no solution. Define 



:= d/^j. Note that is the order of Sj. 



condition for solvability of (|B1|) is or rjyj 
for every k = 1 . 
We have rjSj 



A necessary 
mod d, 



m. We show that this condition holds. 
mod d. From the stabilizer phase 
condition (|15|l . it follows that 



VI. CONCLUSION 



(rj — 1 );•,(/,/)', , + r-j/j = mod 2d, 1 < j < n, 
r jfj — moc ' 2d, j > n, 



We have shown that for the Pauli group, the Clifford 
group and stabilizer states, straightforward extensions in 
Hilbert spaces of arbitrary dimensions can be compactly 
described with matrices over Z^. We have given a way 
of efficiently decomposing an n-qudit Clifford operation 
in 0(n 2 ) one and two-qudit operations. With these tools 
in modular arithmetic, we provide an expansion of an 
arbitrary stabilizer state of n qudits in the standard basis. 



APPENDIX A: THE SMITH NORMAL FORM 

The Smith normal form is a canonical diagonal form 
for equivalence of matrices over a principal ideal ring R. 
In this paper we consider matrices over Id- For any A G 
Z^ xm there exist invertible matrices K G Z ? d ix ™ and L G 



and by definition of y, consequently r^j/j = mod g?, j 
1 ... TO. 

An equivalent system to (|B1|1 is now 



R 



fjj 



mod 



Si 



We define the map b : x = 



[xi 



mod 



b(x) 

T 



which is a homomorphism from the group of vectors of 
length n with entries X4 modulo q^, i = 1 . . . n to the 
group of vectors of length to with entries y'j modulo 
qj/gj, j = 1...TO. IBlfl has a unique solution if b is 
an isomorphism. We prove this by showing that the 
number of elements in both groups are the same and 
that only is in the kernel. It follows from (|14|l and the 
fact that, by definition, the columns of 5 generate a set 
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of d™ elements, that the product of the orders of the 
columns of S is equal to d n , or YYJLi r i ~ d n ■ Therefore, 



TT _ d" 1 ' 

1 I TT m 1 I 



n«=n® 

i=l 



thus the number of elements of both groups are the same. 
Next we show that B T x = mod q if and only if x = 
mod q. We rewrite this as 



Vx G 1 n d : (3v <E : B T x = Q T v mod d) 

(3a;' G Z™ : x = Qa 7 mod d) . 



(B2) 



Proof: 

<=) Q T B is symmetric modulo d. We therefore have 
B T x = B T Qx' = Q T Bx' mod d, so v — Bx' mod d. 

We show that the number of x G satisfying 
the lhs of i|B2(l is equal to the number of x satisfy- 
ing the rhs. The number of elements generated by the 
columns of a matrix is equal to the product of the or- 
ders of the diagonal elements of its Smith normal form. 
Therefore the columns of S T , like the columns of S, 
also generate d n elements. Consequently, the mapping 
s : a G Z^" — > S T a G Z™ is a homomorphism from l? d n 
to a group Y C Z™, with |Y| = tf\ The kernel in Z 2 d n of 
y| = d n elements. Equivalently, with 
[v T w T ], s is a homomorphism from Z^ x ZJJ to Y: 



s contains \Z d n 





V 


) = s T 


V 




w 




w 



Q T B? 



Q T v+B T w. 



There are exactly d n different pairs (v, w) that satisfy 
Q T v + B T w = mod d. Replacing w by —x, we have 
exactly d n pairs (x,v) satisfying B T x — Q T v mod d. 
Fixing such an x, we have a total of Y[ i= i 1i different v 
for which, together with x, the equality still holds (this 
is because we can add an arbitrary multiple of d/qi to 
Vi). Therefore, the total number of x for which a v exists 
such that B T x = Q T v mod d, is equal to <i n /n™=i 
This is equal to the number of x that can be written as 
x = Qx'. □ 
Next, we describe a method for easily finding the solu- 
tion x* of iPJl. We define a diagonal matrix Z G Z™ xm 
with diagonal entries equal to d/q^, k = 1 ...to. I|18() 
is equivalent to the equation ZB T x = Zy mod d. We 
calculate the Smith normal form F = KZB T L mod d. 
Defining ,t' := L~ l x mod d and y' := ifZy mod d, we 
have the following equation Fx' — y' mod d, for which a 
solution x*' £ ZJJ can be easily found (note that this so- 
lution is most likely not unique). We then find x* = Lx* 
mod q. 



APPENDIX C: SIMPLIFICATIONS FOR ODD d 

In this section we consider the special case of odd d. 
Most of the formulas in this paper can be simplified for 



odd d. We will only give an overview and omit the deriva- 
tions, as they are completely analogous to the general 
case. If d is odd, then 2 has an inverse in Z^, equal to 
which we will denote by 2 _1 . 
For odd d, we can use a restricted definition for the 
Pauli group: it contains all d 2n tensor products J2J with 
an additional complex phase factor uo s (instead of a power 
of C). Eq. Q| becomes 



uj s XZ(a)uj € J<Z{b) = lo 



S+e+a T Ub 



XZ(a + b). 



(CI) 



The order of an arbitrary element of this newly defined 
Pauli group is never equal to 2d. In the same way as for 
the general case, we find the image w £ ^(6) of w s XZ{a) 
under conjugation by a Clifford operation, which is now 
defined by C and g = |: 

b = Ca mod d, 

e = S+(g-2- 1 V dmg (C T UC)) T a+ (C2) 
2- 1 a T (C T UC - U)a mod d. 

Note that, contrary to the general case, g is a vector in 
Z^ n . There is no longer a restriction on g. Indeed, from 
, it follows that h in the general setting is always even 
for odd d. Symplecticity of C is of course still required. 
The product of two Clifford operations Q" = Q'Q corre- 
sponds to C" and g" , where 

C" = C'C modd, 

g" = g + C T g' + 2- 1 (v diag (C T (C' T UC'-U)C)- 
C ,T V diag (C" T C/C")) mod d. 

(C3) 

The inverse of a Clifford operation Q defined by C 
and g is defined by C and g' , where 



C = C- 1 = -PC T P mod d, 

9' = -C- T g + 2- 1 (c- T V diag (C T UC)- 

V^giC-TUC- 1 )) modd. 



(C4) 



The definition of a stabilizer state remains the same 
except for the fact that now the stabilizer is a subgroup 
of the restricted Pauli group. Note that for odd d, no 
subgroup of the general Pauli group can be found that 
fulfills all stabilizer conditions but is not a subgroup of 
the restricted Pauli group. Thus, nothing is lost by re- 
stricting the definition of the Pauli group for odd d. A 
generating set for the stabilizer S consists of elements 
u bk XZ{S k ), k = 1 ...to, where S k G Z 2 d n and b k G Z'j 1 . 
Analogously to the definition of g, b is equal to half the 
value of / in the general setting (as it is the exponent of 
lo instead of Q . The stabilizer phase condition Ijlll) on b 
simplifies to: 

Vr G Z™ | Sr = mod d : 
(2b - V dillg (S T US)) T r + r T (S T US)r = mod d. 

(C5) 
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A stabilizer generator matrix change, by applying an in- 
vertible linear transformation R £ Z™ xm to the right on 
S, results in 

S' = SR mod d, 

b> = R T (b-2- 1 V di , e (S T US))+ (C6) 
2- 1 V diag (R T S T USR) mod d. 

A stabilizer state defined by S and b, operated on by a 
Clifford operation defined by C and <?, is a new stabilizer 
state defined by 

S' = CS mod d, 

b' = b + S T {g-2' 1 V diag (C T UC)) + (C7) 
2" 1 V diag (S T (C T UC - U)S) mod d. 

It is not hard to verify that part (ii) of Theorem ^ sim- 
plifies to 

\^) = u; tTMt +P Tt \T(Qt + x*)) (C8) 
tez S 

where M := 2~ 1 QB mod d, 

p := V - V d i a g(M) + B T x* mod i 
In this setting, x* £ Gq :— 7L qi x . . . x Z 9n is defined as the 
unique solution of B T x = —b' mod q. For calculating x* 
we refer to Appendix iBl 
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